NASA— TM— 1 12475 


A/^ & 


A'-o 


/// -3 7 - 1 

/C. v , „r/t- 

Reprinted from 


MATERIALS 
SCIENCE & 

ENGINEERING 

A 



Materials Science and Engineering A208 (1996) 126-130 


Cyclic fatigue of brittle materials with an indentation-induced 

flaw system 

Sung R. Choi a ’\ Jonathan A. Salem b 

a Cleveland State University , Cleveland , OH 44115 , USA 
h NASA Lewis Research Center , Cleveland , OH 44135 , USA 

Received 1 May 1995; in revised form 10 October 1995 



ELSEVIER 


MATERIALS SCIENCE AND ENGINEERING A 

The journal provides an international medium for the publication of theoretical and experimental studies and 
reviews of the properties and behavior of a wide range of materials, related both to their structure and to their 
engineering application. The varied topics comprising materials science and engineering are viewed as appro- 
priate for publication: these include, but are not limited to, the properties and structure of crystalline and 
non-crystalline metals and ceramics, polymers and composite materials. 


Editor-in-Chief 

Professor H. Herman 

Associate Editors 

M, Koiwa (Japan) 

G. Kostorz (Switzerland) 

Editorial Board (MSE A) 

J, Agren ( Sweden ) 

G. Ananthakrishna (India) 

R. J. Arsenault (USA) 

D. Brandon (Israel) 

H. K. D. H. Bhadeshia (UK) 

J. Cadek (Czech Republic ) 

J, B, Cohen (USA) 

J. Driver (France) 

J. D Embury ( Canada ) 

Y. Estrin (Australia) 

H. Fischmeister (Germany) 

C. Garcia de Andres (Spain) 

H. Gleiter (Germany) 

M. W. Grabski (Poland) 

M. Kato (Japan) 

Y. G. Kim (Korea) 

C Laird (USA) 

J. Lendvai (Hungary) 

W. Mader (Germany) 

M. McLean (UK) 


L. Priester (France) 

S. Sampath (USA) 

V. K. Sarin (USA) 

P. Shen ( Taiwan ) 

M. Suery ( France ) 

S. Suresh (USA) 

N. S. Stoloff (USA) 

M. Taya (USA) 

A. K. Vasudevan (USA) 

A. Vevecka (Albania) 

B, Wilshire (UK) 

M. Yamaguchi (Japan) 

T. S. Yen (China) 

Print and Electronic Media Review Editor 

A. H. King (USA) 

Administrative Editor 

Barbara Herman 

Advisory Board (MSE A and B) 

H. Herman, Chairman (USA) 

H. Curien (France) 

M. E. Fine (USA) 

A. Kelly, FRS (UK) 

H. Mughrabi (Germany) 

H. Rangu (Japan) 


Types of contributions 

Original research work not already published, plenary lectures 
and or individual papers given at conferences; reviews of 
specialized topics within the scope of the journal; engineering 
studies; letters to the editor. 

Subscription Information 1996 

Volumes 205-221, each volume containing 2 issues, are 
scheduled for publication. Prices are available from the 
publishers upon request. Subscriptions are accepted on a pre- 
paid basis only. Issues are sent by SAL (Surface Air Lifted) mail 
wherever this service is available. Airmail rates are available 
upon request. Please address all requests regarding orders and 
subscription queries to 

ELSEVIER SCIENCE SA 

P O. Box 564, 1001 Lausanne, Switzerland 

Telephone: (21) 3 20 73 81 

Telex: 450620 ELSA CH Telefax: (21) 3 235 444 

Issues are sent by surface mail after air delivery to Argentina, 
Australia, Brazil, Canada, China, Plong Kong, India, Israel, Japan, 


Malaysia, Mexico, New Zealand, Pakistan, Singapore, South 
Africa, South Korea, Taiwan, Thailand and the USA. Airmail rates 
for other countries are available on request. 

For advertising rates apply to the publishers. A specimen copy 
will be sent on request. 

US and Canadian customers may obtain information from the 
following 

ELSEVIER SCIENCE INC. 

Attn.: Journal Information Center, 655 Avenue of the Americas 
New York, NY 10010, USA. 

Telephone: (212) 633 3750 Telex: 420 643 AEP Ul 
Telefax: (212) 633 3764 

Abstracting and/or Indexing Services 

American Ceramic Society, Cambridge Scientific Abstracts; 
Chemical Abstracts, Current Contents; Engineering Index; FIZ 
Karlsruhe; Fluid Abstracts; Fluidex; Glass Technology Abstracts; 
Inspec Physics Abstracts; Metals Abstracts; Pascal (Centre 
National de la Recherche Scientifique), Physikalische Berichte; 
Research Alert 1 w ; Science Citation Index. 


ELSEVIER 


Materials Science and Kngineering A208 (19%) 126 130 


ENGINEERING 


Cyclic fatigue of brittle materials with an indentation-induced 

flaw system 

Sung R. Choi 1 *, Jonathan A. Salem * 1 

''Cleveland State University, Cleveland , OH 44115, USA 
h NASA Lewis Research Center , Cleveland . OH 44145, USA 

Received 1 May 1995; in revised form 10 October 1995 


Abstract 

The ratio of static to cyclic fatigue life, or 7r ratio’, was obtained numerically for an indentation flaw system subjected to 
sinusoidal loading conditions. Emphasis was placed on developing a simple, quick lifetime prediction tool. The solution for the 
h ratio was compared with experimental static and cyclic fatigue data obtained from as-indented 96 wt.% alumina specimens tested 
in room-temperature distilled water. 

Keywords: Cyclic fatigue; Static fatigue; Indentation-induced flaw; Ceramics 


1. Introduction 

For glass and ceramic materials which have slow crack 
growth (stress corrosion) as the unique, time-dependent 
failure mechanism, it is possible to predict the fatigue life 
under one loading condition from another. Prediction of 
cyclic fatigue lifetime from static lifetime for Griffith flaw 
system can be done analytically by using the ratio of 
static fatigue to cyclic fatigue life, introduced by Evans 
and Fuller [1]. In the case of the indentation flaw system, 
the solution requires a numerical approach since an 
additional driving force appears in the net stress intensity 
factor. This term is attributed to residual contact stress 
produced by the elastic/plastic deformation of indenta- 
tion [2]. Analyses of dynamic and static fatigue of the 
indentation flaw system were carried out previously for 
both postthreshold [3,4] and subthreshold [5] indentation 
flaws. However, no general solution for the cyclic fatigue 
of such a flaw system has been found, although some 
cyclic fatigue data on ceramic specimens containing 
indentation cracks exist [6-9], 


* Corresponding author. Address correspondence to: NASA Senior 
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The purpose of this study is to analyze cyclic fatigue 
of the indentation flaw system under sinusoidal loading 
conditions so that lifetime prediction from one fatigue 
condition to another is readily feasible. For this purpose, 
the complete solution of the ratio of static to cyclic 
fatigue life was obtained numerically in conjunction with 
fatigue parameter (n) and stress ratio (R ratio). The 
solution was compared with experimental data that were 
obtained from static and cyclic fatigue testing of indented 
alumina specimens in room-temperature distilled water. 
The analysis was carried out for material exhibiting a flat 
R curve and a power-law crack propagation as a delayed 
failure mechanism. 

2. Analysis 

In many cases slow crack growth of glass and ceramic 
materials under Mode I loading conditions is described 
by the following empirical, power-law relation [4] 


where i\ a, and t are crack velocity, crack size, and time, 
respectively. A and n are the fatigue parameters which 
depend on material and environment. K x is the 
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Mode I stress intensity factor, and K lc is the Mode I 
critical stress intensity factor or fracture toughness of the 
material with a flat R curve. 


2.1. Natural flaws 


In static fatigue testing a constant applied stress (a) is 
employed. Noting that n > 10 for most glass and ceram- 
ics, one can derive the following static fatigue equation 
[ 10 ] 


tfs = BS'? 


( 2 ) 


where t (s is the time to failure, 5, is the inert strength, and 
B is the parameter associated with fracture toughness, 
crack geometry and fatigue parameters. 

In cyclic fatigue testing a time-varying, periodic stress 
(<x(t)) is applied. The time to failure (f fc ) in cyclic fatigue 
can be expressed as [1] 

/fc = BS'r 2 (3) 

- [Aor dt 
T Jo 

where rx max is the maximum applied stress, /(/) is a 
periodic function defined as a(t) = <r mAX f(t) with a range 
of 0 <f(t) < 1, and r is the period. The ratio of static 
fatigue to cyclic fatigue life, /?, with a condition of o in 
static loading equal to <r max in cyclic loading (<r = <r max ) 
can be obtained from Eqs. (2) and (3) 


A = - = - [V(Ofd/ 

tfc T Jo 


(4) 


Sinusoidal loading is the most common and popular 
wave form used in cyclic fatigue testing. The periodic 
function /(/) for a sinusoidal wave is expressed as 
/(0 = [(1 +/f)/2 + [(l — R)/2]smcot], where/? is the stress 
(or load) ratio, defined as R = a m in /V max with <r min being 
the minimum applied stress, co is the angular velocity. Eq. 
(4) can be solved either analytically or numerically for 
the sinusoidal wave. The h ratio for any other periodic 
loading wave form such as trapezoidal, triangular or 
square can be solved analytically with a straightforward 
procedure [1,11]. Eq. (4) was solved analytically for a 
sinusoidal wave by using the series expansion technique 
of Evans and Fuller [1]. Their analytical solution is 
expressed as 




R 

2(1 +~R)_ 


1 +R\” 


(5) 


Eq. (5) is still complicated to solve and valid only for 
integer values of n. Therefore, it is desirable to solve Eq. 
(4) numerically to cover any integer or real value of n for 
a full range of R = 0- 1 .0, even though some limited data 
on the h ratio exist [1,12]. 

Fig. 1 shows the results of the numerical solution of 
Eq. (4) as a function of n for R ratios from R = 0 to 
R= 1.0. If log h is treated as a linear function of log n , 
coefficients of correlation of r coef > 0.9966 result. Thus it 



FATIGUE PARAMETER, n 

Fig. 1. Numerical solution of the h ratio as a function of fatigue 
parameter (n) for different levels of R ratios with the natural flaw 
system. Each line represents a best-fit regression line. 

is possible to obtain an approximate (but accurate) 
relation between log h and log n for a given R ratio by 
using the regression analysis 

log h = a log n - h /? (6) 

where a and p are the regression coefficients. 

Table 1 shows the regression coefficients a and p. For 
n > 10, the maximum error in h ratio, as compared with 
the exact solution, was 1.0% for R = 0 to 0.7, and 4.5 % 
for R = 0.8 to 0.9. Therefore, the table gives a convenient 
and accurate means of determining the h ratio for a 
natural flaw system for any given value of n , either integer 
or real. The complication of a lack of a simple analytical 
solution to Eq. (5) for real numbers can be eliminated 
with Eq. (6) and Table 1. 

2.2. Indentation flaws 

Because of the residual contact stress produced by 
elastic/plastic indentation deformation [2], an addi- 


Table 1 

Regression coefficients of the h ratio (Eq. (6)) for natural and 
indentation flaw systems 


R 

Natural flaws 


Indentation flaws 

0L 

a 

0L 

p 

0.0 

-0.4939 

-0.2607 

— 

— 

0.1 

-0.4955 

-0.2348 

-0.4856 

-0.1925 

0.2 

-0.4976 

-0.2051 

-0.4885 

-0.1612 

0.3 

-0.5005 

-0.1703 

-0.4925 

-0.1240 

0.4 

-0.5049 

-0.1282 

-0.4977 

-0.0804 

0.5 

-0.5112 

-0.0765 

-0.5049 

-0.0268 

0.6 

-0.5195 

-0.0115 

-0.5126 

-0.0383 

0.7 

-0.5291 

0.0705 

-0.5198 

0.1172 

0.8 

-0.5353 

0.1743 

-0.5173 

0.2075 

0.9 

-0.5138 

0.2977 

-0.4743 

0.2932 

1.0 

0.0000 

0.0000 

0.0000 

0.0000 
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tional term appears in the net stress intensity factor and 
an analytical solution of cyclic fatigue for the indenta- 
tion flaw system is not feasible. The solution needs to 
be done via numerical methods, as done previously for 
both dynamic and static fatigue analyses involving in- 
dentation fracture mechanics [3,4], As in the previous 
studies [3,5], normalized variables are introduced here 
as follows: 


K* = 


K x 


J = — t ; <t* 

ct m 


C* = — 
a m 


(7) 


where K *, 7, a* and C* are, respectively, normalized 
stress intensity factor, normalized time, normalized 
maximum applied stress and normalized crack size. <r m 
and a m are, respectively, the strength and critical crack 
size in the inert condition. Using these variables, the 
crack growth rate of Eq. (1) and net stress intensity 
factor are expressed as follows: 


dC* 

-d7 = [ ** ] 


n 


= - + ^ C*- 32 


(<- 


... ri + R i -r . t 

<ri = [— + — s - n ( 


( 8 ) 


Note that the normalized net stress intensity factor 
consists of two terms: the first is a function of the 
remote applied stress and the second is related to the 
residual contact stress [2,3]. The differential form of Eq. 
(8) was solved numerically using a fourth-order 
Runge-Kutta method. The initial condition was C* = 
0.3967 at 7 = 0 and the instability conditions were 
K* = 1 and dtf */dC* > 0. The solution was initiated to 
determine the normalized time to failure (7 f ) as a func- 
tion of normalized maximum applied stress for the 
selected values of n = 5-160. This procedure was con- 
tinued for the range of stress ratios from R = 0.1 to 1.0. 
Since the solution is independent of frequency (/) as 
long as (oaj A > 2nf, any particular values of coaJA 
satisfying coaJA > In can be chosen. A value of coa m / 
A = 100 was used in this analysis. 

Fig. 2 shows a typical example of the normalized 
time to failure as a function of normalized maximum 
applied stress for stress ratios of R = 1.0, 0.5 and 0.1. 
Fatigue susceptibility increases with increasing R ratio, 
yielding a maximum at R= 1.0 (static fatigue). Note 
that regardless of R ratio the curves converge to a * = 
1.0, in which the inert strength with no slow crack 
growth is defined. Similar to the previous static fatigue 
analysis [4], the slope in Fig. 2 is not representative of 
a 'true’ fatigue parameter of /?, due to the effect of 
residual contact stress. The slope in Fig. 2, which is 
called ‘apparent’ fatigue parameter ( n' ), was found to 
have the following approximate relation 


n = 4/3n' — 2/3 


(9) 



Fig. 2. Typical example of numerically obtained normalized failure 
time (7 f ) as a function of normalized maximum applied stress (<r*) for 
three different R ratios of /? = 0.1, 0.5 and 1.0 for the indentation 
flaw system. 


which reduces to the relationship in static fatigue of 
indentation cracks [4]. This indicates that the relation- 
ship between the true (n) and apparent (n f ) fatigue 
parameters is independent of R ratio. 

Based on the results as shown in Fig. 2, the h ratio 
was determined for a given R using the relation h = JJ 
7 fc ( = hs/hcA where 7 fs and 7 fc are the normalized time 
to failure corresponding to static and cyclic fatigue, 
respectively. Fig. 3 shows a summary of the h ratio as 
a function of n for R ratios from /? = 0.1 to R= 1.0, 
where log h was plotted against log n as for the natural 
flaw system (Fig. 1). Again the excellent coefficients of 
correlation of r coef >0.9911 allow one to obtain a rela- 
tionship based on Eq. (6). The results of such regression 
analysis for a and are shown in Table 1. The maxi- 



Fig. 3. Numerical solution of the h ratio as a function of fatigue 
parameter ( n ) for different R ratios with the indentation flaw system. 
Each line represents a best-fit regression line. 
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mum error associated in the h ratio was observed to be 
about 4%, occurring at R = 0.9 for n > 10. Hence, 
lifetime prediction of an indentation flaw system with one 
fatigue loading condition can be done quickly and 
accurately from another by using Eq. (6) in conjunction 
with Table 1. 


2.3. Evaluation of fatigue parameters for indentation 
flaws 


The “true” fatigue parameter n for both static and 
cyclic fatigue can be determined using Eq. (9) once the 
“apparent” fatigue parameter n is determined from the 
slope of fatigue data. The fatigue parameter A in static 
fatigue was determined previously [4] 


A = 



( 10 ) 


where /' is the intercept of static fatigue data, expressed 
as t ts <r n = a' s . Since the following relation for the inden- 
tation flaw system holds 


A 

A 


A s (7 




( 11 ) 


where a c is the intercept of cyclic fatigue data as 
ffcfl’max = K- Hence, from Eqs. (10) and (11), the parame- 
ter A in cyclic fatigue is obtained 

f2n\ 2 .1 

A= — (T"a m — (12) 


3. Experimental procedure 

Static and cyclic fatigue tests of indented alumina (96 
wt.%, ALSIMAG 614, General Electrical Ceramics) 
flexure specimens were carried out in room-temperature 
distilled water using a four-point bend fixture of 6.05 mm 
inner and 19.05 mm outer spans. The nominal dimen- 
sions of the test specimens were 4 mm by 5 mm by 25 
mm, respectively, in height, width, and length. The center 
(5 mm side) of each specimen was indented in air for 
about 20 s using a Vickers microhardness indenter 
(Zwick, model 3212, Germany) with an indentation load 
of 49 N. Static fatigue testing was conducted by loading 
indented specimens in a lever-arm creep machine (ATS) 
with constant stress levels of 95 to 133 MPa. Cyclic 
fatigue testing of indented specimens was conducted by 
sinusoidal loading with a stress ratio of R = 0.5 and a 
frequency of /= 5 Hz with a servohydraulic testing 
machine (Instron, Model 8562). The maximum applied 
stress in cyclic fatigue ranged from 100 to 130 MPa. The 
sinusoidal wave shape was frequently checked and ver- 
ified by a digital storage oscilloscope. 

4. Results and discussion 

The results of the static and cyclic fatigue testing of 
the indented 96 wt.% alumina specimens in room-tem- 



(MAX) APPLIED STRESS, o [MPa] 

Fig. 4. Results of static and cyclic fatigue testing obtained from 
indented 96 wt.% alumina flexure specimens tested in room-tempera- 
ture distilled water. The solid lines represent the best-fit lines in a plot 
of log r, vs. log <r mav 

perature distilled water are depicted in Fig. 4. Consis- 
tent with the results shown in Fig. 3, the alumina is 
slightly more susceptible to fatigue in static loading 
than in cyclic loading within the experimental range 
used. The “apparent” fatigue parameter was obtained 
from the data as ri = 32.92 ± 6.01 and 35.37 ±4.79, 
respectively, for static and cyclic fatigue. The corre- 
sponding “true” fatigue parameter was obtained from 
Eq. (9) as n = 43.23 ± 8.01 and 46.49 ± 6.39, respec- 
tively, for static and cyclic fatigue. The prediction of 
cyclic fatigue from static fatigue data can be done using 
Eq. (6) together with Table 1. The prediction thus made 
is presented in Fig. 5, where the cyclic fatigue data 
obtained from the experiments are included for com- 
parison. The prediction somewhat overestimates (less 


w 

as 

tD 

-j 

< 


o 


w 




e- 
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Fig. 5. Prediction of cyclic fatigue lifetime from static fatigue data for 
indented 96 wt.% alumina specimens tested in room-temperature 
distilled water. 
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than one order of magnitude) the actual cyclic time to 
failure. 

It has been reported that, for certain ceramic materi- 
als, damage accumulation and/or fatigue synergisms 
can be active in cyclic fatigue, resulting in more fatigue 
susceptibility in cyclic than in static loading [8,13-16]. 
If such synergies exist or if creep at elevated tempera- 
tures is simultaneous with cyclic fatigue, the analysis 
given in this paper may not be valid. However, the 
analysis may still provide, by comparing with experi- 
mental data and by using fractographic analysis, some 
clues with which the prevailing mechanism associated 
with failure can be pinpointed. It should be noted that 
most of glass and ceramic materials under static fatigue 
loading conditions in a room-temperature moisture en- 
vironment are subjected to one failure mechanism, 
stress corrosion. Therefore, in view of the reasonable 
agreement between the static and cyclic fatigue data 
shown here, it can be stated that stress corrosion is a 
governing failure mechanism, either in static or in cyclic 
loading conditions. This is supported by the fact that 
the material contains a large amount of glassy phase, 
which is highly susceptible to stress corrosion in an 
aqueous environment. The fatigue parameters for this 
material system, therefore, can be obtained by using 
either static or cyclic fatigue testing. In terms of testing 
economy, however, static fatigue is preferred because of 
the much higher testing costs associated with cyclic 
fatigue testing. 

As mentioned before, the analytical solution of the h 
ratio for other loading functions such as trapezoidal, 
triangular and square wave forms can be easily ob- 
tained for the natural flaw system because of their 
mathematical simplicity [1,11]. However, the conven- 
tional numerical solution for these loading functions 
may not be simple for the indentation flaw system, 
since the functions, unlike the sinusoidal wave, are not 
continuous but discontinuous in nature. In this case, 
each discontinuous periodic function should be con- 
verted into a continuous periodic function so that nu- 
merical solution is feasible. This may be done by using 
the Fourier analysis from which a discontinuous func- 
tion can be approximated into a continuous function. 


4. Conclusions 

The ratio of static to cyclic fatigue life, or h ratio, 
was obtained numerically with an emphasis on the 
indentation flaw system subjected to sinusoidal loading 
conditions. The h ratio decreases with increasing n and 
decreasing R ratio. The solution provides a simple and 
quick methodology to predict the lifetime of one fatigue 
condition from another for indentation flaw systems (as 
well as the natural flaw system) when the governing 
failure mechanism is stress corrosion. The solution was 


compared with static and cyclic fatigue data obtained 
from indented 96 wt.% alumina specimens tested in 
room-temperature distilled water. 
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